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• • ' 1 INTRODUCTION 



ABSTRACT 

We have analysed archival data taken with the Australia Telescope Compact Array 
(ATCA) during 2001-2003 and detected nine new interstellar and circumstellar H2O 
masers in the LMC. This takes the total number of star formation H2O masers in the 
LMC to 23, spread over 14 different star forming regions and three evolved stars. Three 
H 2 maser sources (N105a/MC23, N113/MC24, N157a/MC74) have been detected in 
all the previous observations that targeted these sites, although all show significant 
variability on timescales of decades. The total number of independent H2O maser 
sources now known in the LMC means that through very long baseline interferometry 
astrometric measurements it will be possible to construct a more precise model of the 
galactic rotation of the LMC and its orbital motion around the Milky Way Galaxy. 

Key words: masers - starsdormation, mass-loss - Magellanic Clouds. 
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There are more than a thousand H2O maser sources 
currently known in the Milky Way (M W) galaxy (e.g., 
IValdettaro et al. Il200il ; I Walsh et al. 112011! ). In contrast, the 
total number of maser sources detected in all species in the 
Large and Small Magellanic Clouds (LMC, SMC) is around 
20 fsee lEllingen et al. II20T0I , hereafter EBCQF10, and refer- 
ences therein) . The maser sources in the LMC and the SMC 
are important objects for addressing a range of astrophysical 
questions. 

H2O masers are typically associated with outflow activ- 
ity from young stellar objects (YSOs) and the final copious 
stellar mass loss phase of red giant and supergiant stars. All 
of the H2O masers found in the Magellanic Clouds (MCs) 
to date are associated with either massive YSOs or red su- 
pergiants, both of which trace the present sites of star for- 
mation in these galaxies. Taking into account the relatively 
low metallicity of the MCs compared to the MW, some dif- 
ferences are expected in the dominant physical mechanisms 
and the timescales that govern the process of star formation 
and stellar mass loss. 



Very long baseline interferometric (VLBI) observations 
can resolve each maser source into a cluster (or clusters) of 
compact maser features, whose three-dimensional motions 
(line-of-sight velocities and proper motions) can be mea- 
sured. With current instruments with which astrometric ac- 
curacy of 10 micro- arcseconds (/Lias) is achievable, such stud- 
ies are possible for star forming r egions in nearby galaxie s 
at distances of up to 1 Mpc (e.g. iBrunthaler et alT1l2005l ). 
Measuring the motions of the maser sources with respect to 
nearby (in terms of angular separation) quasars, it is possi- 
ble to determine both a space motion and a trigonometric 
parallax distance (given sufficient astrometric accuracy) of 
the maser source. Within the MW on kilo-parsec scales, H2O 
maser sources are confined to the MW thin disk and have 
been used to derive fundamental Galactic s tructure param- 
eters jReid et al. 1120091 ; iHonma et al. Il2012[ ). The LMC and 
SMC are at dis tances of approxim ately 50 and 60 kpc re- 
spectively fe.g.. lCioni et ain feOOO). so the amplitude of the 
trigonometric parallax is approximately 20 pas. This is too 
small to enable accurate parallax distance det erminations 
with current instruments and techniques (e.g. iReid et al. I 
2009), but may be feasible through statistical approaches, 
or with future space-based VLBI missions. 



* E-mail: hiroimai@sci.kagoshima-u.ac.jp 



With current, ground-based VLBI we are able to deter- 
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Table 1. Setup status of the analysed data of the ATCA observations towards the LMC. 



Program 


Date 






Bandwidth 






code 


(yyyy/mm/dd) 




Arrayf 


Scan timet (MHz) 


Aspect § 


Targets 


C901 .. 


2001/01/07 


3 


750C 


18-19 X 3-4 8 


1024 


N105a, N113, N157a, N159, N160a, and 0506-612 


C973 .. 


2002/11/01-07 


3 


Wnnn^[ 


1-3 x 2-3 16 


513 


186 sources in HII regions and O-rich AGB stars 




2003/05/01-04 


6 


EW352 


5-7 x 2-3 16 


513 


28 sources in HII regions and O-AGB 815 



*Number of scans X duration per scan in min. t Array configuration. ^Number of available baselines. §Number of spectral channels. 
KThe array configurations are named as combination of antenna positions labelled as W plus three digit numbers such as: W100, 
WHO, W147, W168 and W196. 



mine the space motion of maser sources within the MCs by 
measuring and correcting for the intrinsic internal motions 
of maser features in the source. The locations and past orbits 
of the MCs are important and controversial topics in stud- 
ies of the star formation history and process of interactions 
of th e galaxies in the MW-MC system fe.g.. ll5iaz fe Bekkil 
120121 ). Astrometric observations using data from ground- 
based optical telescopes and the Hubble Space Telescope 
have been used to measure the proper motions of the MCs 
JPiatek et al. 1 120081 : IVieira et al~l l20ld : iKallivavalil" et al. 
120131 ). However, there exist non- negligible discrepancies 
among the results. The optical astrometry results are af- 
fected by the inclusion of stars from a variety of popula- 
tions with different dynamical characteristics located along 
the same lines of sight. They are also affected by the galac- 
tic ro tation model adopted (e.g.. lHaschke. Grebel fc Duffaul 
I2012T I, which is required to extract the centre-of-mass space 
motions of the galaxies from the observed proper motions. 

In order to analyse the 3-D kinematics of the individual 
H2O maser sources in the MCs and to derive the dynamical 
parameters of these galaxies, efforts to increase the num- 
ber of identified H2O masers in these galaxies are impor- 
tant. Statistical analysis of the 3-D kinematics of the maser 
sources will enable us to compare the kinematic properties 
of the MCs with those of the MW. In order to construct 
dynamical models of the MCs, each galactic model requires 
the present centre-of-mass space motion, the velocity field of 
the galactic rotation and the distance to the galaxy as free 
parameters in the model fitting. A large number of maser 
sources (>20) are required to make such an analysis feasible. 
In this paper, we present H2O masers in the LMC detected 
with the Australia Telescope Compact Array (ATCA), in- 
cluding nine newly discovered since the work of EBCQF10. 



useable. Data reduction was undertaken using the MIRIAD 
package, following the standard procedures for ATCA data. 
The visibility amplitude and phases were calibrated by ref- 
erencing them to observations of flux density and phase cal- 
ibrator sources. Phase calibration solutions were obtained 
from the scans on the phase calibrators, J045005. 4— 810102 
and J050644. 0—610941, one of which was observed every 10- 
20 min. Although there was a difference in the channel reso- 
lution for the correlator outputs in the different sessions, all 
data were smoothed to the achieve the same velocity resolu- 
tion of ~0.5 km s _1 . We then constructed synthesis images 
for each source to search for H2O masers and to check the 
image fidelity. 

For sources where a maser detection was confirmed and 
there was sufficient coverage of the uv-plane (mainly in the 
C973 2003- June session), we attempted to obtain the spec- 
trum and the coordinates of the maser emission from the 
image cube. Even for the sources with sufficient uv-coverage, 
there existed significant side lobes, making the astrometry 
difficult. The angular resolution of the ATCA observations 
was typically ~0".5 (corresponding to a linear resolution of 
~0.1 pc at the distance of the LMC), within which most ve- 
locity components of a single maser source will be contained. 
Because we are mainly interested in the number of indepen- 
dent maser sources rather than the internal structures of the 
individual maser sources in this paper, we have only deter- 
mined the coordinates of the brightest velocity components 
of the maser sources. The uncertainty in the measured maser 
positions is 1-2" in the cases where we could not uniquely 
determine the brightest point of the maser source due to 
the high side lobe levels. Columns 2 and 3 of Table [5] give 
the coordinates of the detected maser sources where we were 
able to determine them. 



2 ANALYSIS OF ATCA ARCHIVAL DATA 

We analysed archival observations of H2O maser sources in 
the MCs taken in ATCA observations made on 2001 January 
7 (program C901), and in the period of 2002 February-2003 
May (C973fl- Tabled] summarises the observing strategy 
and correlator setup used for the observations. In each of 
the observing sessions all six 22-m antennas were used in the 
array, but for some observations only three baselines were 

1 The C973 observations also targeted sources in the SMC, which 
we also reduced. However, the SMC data have b een independently 
analysed and published by iBreen et al. I ||2013T ) and in this paper 
we focus only on the sources in the LMC. 



3 RESULTS 

A total of twelve 22-GHz H2O masers were detected to- 
wards the LMC in the C901 and C973 ATCA observations, 
eight and one of which are newly discovered interstellar and 
circumstellar sources, respectively. Table [5] lists the sixteen 
interstellar H2O masers in the LMC known to date (see also 
th e latest review of maser source surveys towards the MCs 
in Ivan Loonll2012l ). Table |3] gives the parameters of the cir- 
cumstellar maser newly detected. Figure [1] shows the ATCA 
cross-power spectra of the H2O masers in the LMC. 

The stellar source O-AGB 815 (2MASS 
J05351409-6743558) corresponds to an M4-type vari- 
able star (HV 1001) with a A'— magnitude of 8.14. Its 
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Table 2. Interstellar 22-GHz H2O masers in the LMC. Some of the source names in the first column are those regi ste red in the data 
archiv e. The source coordinates are cited from the latest result f rom t he p resent work, EBCQFl O. lLazendic et al~ 1 l|2002h or fOliveira et al. I 
ll2006h . References for the coordinates are, l: lLazendic et al. I j2002l) ; 2: lQliveira et al. I ll2006h ; 3: EBCQF10; 4: Present work. 1-tr noise 
level s on the images (inclu ding side lobes) are given in the present results. References and catalogs for the corresponding objects are, 
GC: lGruendl fc Chul hoO$ ); 2MASS: Two-micron All Sky Survey Point Source Catalog; 2MASX: 2MASS Extended Source Catalog. 



RA (J2000) Dec. (J2000) Velocity Peak 1-cr 

Range Peak flux noise 



Source name 


h 


m 




S 


/ 




// 


( km s 


- 1 ) 




(Jy) 


Ref. 


(Jy) 


Counterpart stellar objects 


IRAS F04521-6928 . . . 


4 


51 


53 


.90 


-69 23 


27. 


.8 


205- 


-245 


233. 


.4 


2.5 


4 


0.08 


GC J045153. 29-692328. 6 


HII 1107 


4 


52 


09 


.74 


-66 55 


21. 


.8 


250- 


-275 


273. 


.5 


0.9 


4 


0.07 


2MASS J04520916-6655223 


HII 1186 


4 


54 


00 


.10 


-69 11 


54. 


5 


245- 


-255 


245. 


.5 


5.0 


4 


0.08 


2MASS J04540020-6911557 


IRAS 05011-6815 


5 


01 


01 


.83 


-68 10 


28. 


2 


263- 


-266 


264 


A 


0.08 


3 




GC J050101. 80-681028. 5 


N105a/MC23 


5 


09 


52 


.00 


-68 53 


28. 


6 


252- 


-260 


254. 


.6 


1.4 


1, 3 




GC J050952. 26-685327.3 




5 


09 


52 


.74 


-68 53 


21. 


.9 


235- 


-260 


244. 


.2 


1.1 


2, 4 


0.09 


GC J050952. 26-685327.3 




5 


09 


52 


.55 


-68 53 


28. 


.0 


246- 


-268 


266 


.2 


0.8 


3 




GC J050952. 26-685327.3 


N113/MC24 


5 


13 


17 


.17 


-69 22 


21. 


6 


249- 


-286 


257. 


.1 


2.0 


3 




GC J051317.69-692225.0 




5 


13 


17 


.69 


-69 22 


27. 





245- 


-250 


249. 


.1 


0.3 


1, 3 




GC J051317.69-692225.0 




5 


13 


25 


.31 


-69 22 


42. 


.4 


215- 


-257 


237. 


.8 


52.3 


2, 4 


0.11 


2MASS J05132499-6922442 




5 


1.3 


25 


.15 


-69 22 


45. 


1 


228- 


-264 


253 


.1 


116 


1, 2, 3 




2MASS J05132499-6922442 


N113a (in N113) 


5 


13 


22 


.36 


-69 22 


38. 


6 


235- 


-265 


238 


.1 


3.1 


4 


2.56 


GC J051321.43-692241.5 


IRAS 05202-6655 


5 


20 


16 


.90 


-66 52 


53. 


3 


285- 


-300 


287. 


.0 


2.0 


4 


0.06 


2MASX J05201653-6652544 


NGC 1984 


5 


27 


11 


.32 


-69 08 


05 


.9 


255- 


-265 


258 


.3 


.3.3 


4 


0.07 


2MASS J05274010-6908044 


N157a/MC74 


5 


38 


46 


.5.3 


-69 04 


45. 


.3 


253- 


-271 


269 


.0 


6.5 


1, 3 




GC J053848.42-690441.6 




5 


38 


49 


.25 


-69 04 


42. 


.1 


263- 


-268 


267 


.3 


0.5 


3 




GC J053848.42-690441.6 




5 


38 


16 


.80 


-69 04 


45. 


.5 


235- 


-255 


249 


.9 


0.6 


2, 4 


0.09 


GC J053848.42-690441.6 




5 


38 


52 


.67 


-69 04 


37 


.8 


265- 


-269 


266 


.9 


1.4 


3 




2MASS J05385273-6904374 


N160a 


5 


39 


13 


.83 


-69 38 


33 


.8 


236- 


-261 


251 


.6 


1.8 


1, 2, 3 




GC J053943.82-693834.0 




5 


39 


38 


.97 


-69 39 


10 


.8 


252- 


-260 


253. 


.1 


2.5 


3 




GC J053939.02-693911.4 


N159 


5 


39 


29 


.23 


-69 47 


18. 


.9 


242- 


-255 


248 


.3 


1.1 


1, 3 




GC J053929.21-694719.0 


N214b 


5 


39 


53 


.51 


-71 09 


51. 


6 


225- 


-230 


228 


.7 


7.8 


4 


0.09 


2MASS J05395357-7109533 


N180 


5 


48 


26 


.44 


-70 08 


48. 


.8 


130- 


-140 


139 


.8 


0.6 


4 


0.07 


2MASS J05482622-7008502 



Table 3. Same as Table [2] but for a circumstcllar 22-GHz H2O maser in the LMC, which was newly detected in the present work. 



RA (J2000) Dec. (J2000) Velocity Peak 1-cr 

Range Peak flux noise 

Source name h m s o 1 11 ^ ^ m s -i) (jy) (jy) Counterpart stellar object 

O-AGB 815 5 35 14.57 -67 43 54.5 230-236 235.5 0.18 0.08 HV 1001, 2MASS J05351409-6743558 

Note. The coordinates are cited from the Harvard Variable (HV) Star Catalogue. 



H2O maser is apparently weak (~0.2 Jy) and the maser 
luminosity is comparable to or lower than that seen in 
supergiants in the MW (e.g., VY CMa with an H2O maser 
flux density of u p to 1000 Jy at a distance of 1.1 kpc, 
IChoi et al. 1120081 ). 

The eight newly detected interstellar H2O masers 
(those where Column 7 of Table [2] contains only 
the reference 4) are IRAS F04521-6928, HII 1107, 
HII 1186, N113a (in N113/MC24), IRAS 05202-6655, 
NGC 1984/OH 279.6-32.6, N214b and N180. In some re- 
gions, H2O maser sources were spread over an area of order 
1'. Taking into account the position uncertainty, we regarded 
maser emission located within a radius of 1"— 2" as a single 
source and they are shown in a single line in Tabled There 
are fifteen independent H2O maser sources that are sepa- 
rated from each other by more than 1' (15 pc at the distance 
of the LMC). 

Some bright H2O m aser sources were first det ected more 
than 25 years ago (e.g., IScalise Jr. fc Braj|l98lh . However, 
even for regions where there has been persistent maser emis- 
sion over a long period, the sources have exhibited large 



flux density variability and sometimes drop below the detec- 
tion threshold. For example, in the brightest maser source 
N113/MC24, the peak velocity of the brightest maser com- 
ponent has dramatically changed and the strongest spectral 
feature observed in this paper (observation in 2003 May) 
was not detected in the observations of EBCQF10 (2008 Au- 
gust). As a result, apparent position shifts (>l" or >0.25 pc) 
may occur due to the maser sources with such flux variabil- 
ity in the different locations (a Christmas tree effect). Clus- 
tering of H2O maser sources on parsec scales is observed 
in some Galactic star formation regions (e.g., W51 North 
and Main; W3(OH) and W3 IRS5). Therefore, the observed 
position shifts may be due to variation in the relative flux 
density between maser sources within a cluster, although we 
cannot rule out the possibility that these shifts occur just 
due to the relative positional errors of the observations made 
at different epochs. N113/MC24 and N113a present one ex- 
ample of a cluster of H2O maser sources (|Lazendic et al. I 
l2002i ; IOliveira et al. II2OO6] ; EBCQF10). 

We looked for corresponding objects within 10" of the 
maser sources and these are listed in Column 9 of Table 
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Figure 1. H2O maser spectra taken with the ATCA in the 2003 May session, except for N113a taken in the 2002 November session. 



[2] Only the nearest sources from the maser position are 
listed. The corresponding 2MASS sources are located within 
3" of the masers. Note that IRAS 05202-6655 seems to 
be associated with an extended infrared source (2MASX 
J05201653— 6652544) and its population i s unc lear. Prox- 
imity to sources listed bv lGruendl fc C hu (2009) may be a 
good indication of association with a YSO. 



4 DISCUSSION 

Here we estimate possible detections of H2O masers in the 
LMC using an H2O m aser luminosity function (H2O-LF). 
iGreenhill et"ai~l (|l990r i derived the H2O-LF for the MW, 
which will be improved by further unbiased H2O maser sur- 
veys in the MW (e.g.. IWalsh et al. Il20ll ). Similar to the 
approach used bv lDarlind ( 2011 ) to estimate the expected 
number of H2O masers in M31, we need to rescale the H2O- 
LF using the relative star formation rates of t he MW and 
the LMC (~4 and ~ OAM ^yr' 1 , respectively. iDiehl et al. I 
120061 ; ISkibba et al71l2012h . 

The shortest on-source integration time for the present 
data was ~3 min. With a velocity resolution of 0.5 km s _1 , 
the corresponding detection limit for the peak of any H2O 
maser emission is estimated to be 1 Jy. Adopting a distance 
to the LMC o f 51 kpc (corresponding to a distance modulus 
of 18.55, e.g.. lCioni et al. Il2000l ; iHaschke, Grebel fc Duffaul 



our detection limit corresponds to an isotropic H2O 
maser luminosity of 2.4 x 10 -5 Lq. We calculate the ex- 
pected number of H2O masers with the luminosities higher 
than this value in the LMC to be ~15. This is in good agree- 
ment with the actual number of the identified sources over 
1 Jy (16 sources in Table [2]). This supports the argument 
that the number of H2O maser detections in the LMC is 
consistent with MW expectations when the difference in the 
star formation rate s between the LM C and the MW are 
taken into account (|Green et al. I l2008l ; EBCQF10, and ref- 
erences therein). Although the same approach when applied 
to M31 results in an underestimation of the expected number 
of detections (prediction of ~3 compared to the 5 detected, 
iDarlindlaniT ). the difference is less than a factor of 2 and 
within the estimated uncertainty of this approach. 

EBCQF10 used the LMC YSO catalogue of 
iGruendl fc Chul l|2009h to investigate the infrared properties 
and spectral energy distribution (SED) of star formation 
regions with an associated H2O maser and compared them 
to the entire sample. They showed that EhO-associated 
YSOs have higher central mass, larger envelope radius, 
higher ambient density and higher total luminosity than 
the general YSO population. EBCQF10 classified YSO 
within 2"of the H2O maser position as being associated. 
Applying the same criterion to the 8 new LMC H2O 
masers listed in Table H we find three (HII 1186, N214b 
& N180) have a definite YSO each within 2" (table 9 
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Figure 2. Distribution of all interstellar H2O masers identified to 
date in the LMC on a gnomonic projection map. The background 
grey image is the optical image of the LMC taken from Sky View: 
http://skyview.gsfc.nasa.gov 



5 CONCLUSION 

With the current sample of H2O masers in LMC (19 in to- 
tal) and the existence of position-reference quasars within a 
few degrees of arc of these masers (Imai et al. in prepara- 
tion), VLBI astrometric observations to estimate the centre- 
of-mass space motion vector for the LMC are feasible. They 
may provide an opportunity for independently checking the 
results of optical astrometry measurements. Furthermore, 
measurement of the trigonometric parallax of the LMC 
(n ~ 20/ias ) may be possible. With a position measurement 
accuracy achieved by the state-of-the-art VLBI astrometry 
mentioned previously, an individual measurement of the par- 
allax will yield only a marginal detection. However, with 
independent measurements from ~10 regions, each with po- 
tentially a few detections of the maser spot parallax, it will 
be possible to statistically improve the reliability and pre- 
cision of the parallax detection. In the same observations 
it will also be possible to investigate the spatio-kinematical 
structure for each H2O maser source and to study the phys- 
ical properties of the star formation regions and the final 
stellar mass loss in the MCs. 



of iGruendl fe Chi 120091 ). In contrast, for the 16 LMC 
interstellar masers investigated by EBC QF10, 11 were 
within 2" of a definite YSO identified by IGruendl fe Chul 
(2009). The cause of this discrepancy is not clear, it may 
indicate that the newly detected LMC H2O maser sources 
are primarily associated with more evolved star formation 
regions than the sample of EBCQF10, which targeted all 
the CH3OH, OH and H2O masers known at that time. The 
three associations from the newly detected masers have 
similar SED properties to the maser-associated sources of 
EBCQF10. However, none of them are sources identified by 
EBCQF10 as a high-probability maser candidate, so future 
searches towards these sources may provide additional H2O 
masers in the LMC. 

Including the three circumstellar H2O masers associ- 
ated with red supergia nts (IRAS 04553-6825, 05280-6910, 
Ivan Loon et al. Il200ll ; HV 1001) as members of the young 
stellar population, there are 19 independent maser sources 
in the LMC. Even limited to 16 interstellar masers brighter 
than 0.5 Jy, through measurement of the proper motions of 
these masers (better than 100 ^tas yr _1 in a time baseline of 
one year) it is possible to obtain 48 observables (the three- 
dimensional velocity vector components for each source). A 
basic dynamical model of the LMC requires nine free param- 
eters in the fitting: a source distance, a three-dimensional 
centre-of-mass space motion vector, the galactic rotation 
axis inclination and position angle and a third (or higher 
order) polynomial galactic rotation curve. Fig. [2] shows that 
the H2O masers are distributed throughout the LMC, miti- 
gating the likelihood of bias in the model fitting which may 
occur if only a small number of independent locations were 
available. If the distance to the LMC and the axis of the 
galactic rotation are fixed in the model fitting (using previ- 
ous accurate determinations), it will be possible to extract 
the velocity field of the LMC and more precisely determine 
the centre-of-mass space motion of the LMC from maser 
proper motion observations. 
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